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In this paper we present a single-soliton two-component cellular au- 
tomata (CA) model of waves as mobile self-localizations, also known 
as: particles, waves, or gliders; and its version with memory. The 
model is based on coding sets of strings where each chain represents a 
unique mobile self-localization. We will discuss briefly the original soli- 
ton models in CA proposed with filter automata, followed by solutions 
in elementary CA (ECA) domain with the famous universal ECA Rule 
110, and reporting a number of new solitonic collisions in ECA Rule 
54- A mobile self-localization in this study is equivalent a single soliton 
because the collisions of these mobile self-localizations studied in this 
paper satisfies the property of solitonic collisions. We also present a 
specific ECA with memory (ECAM), the ECAM Rule (f>mmaj-A, that 
displays single-soliton solutions from any initial codification (including 
random initial conditions) for a kind of mobile self-localization because 
such automaton is able to adjust any initial condition to soliton struc- 
tures. 

Keywords: soliton, elementary cellular automata, memory, localiza- 
tions, collisions, computability 
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| 1. Introduction 

A soliton can be denned in an informal as follows: when two soli- 
tary waves travel in opposite directions and collide, they emerge after 
collision with the same shape and velocity asymptotically. The phe- 
nomenon of solitary wave first recognized by English engineer John 
Scott Russell [59 ]Jj and first formalized by Diederick J. Korteweg and 
Gustav de Vries in 1895 25J. However, in 1965 the physician Martin 
Kruskal coined the phenomenon of solitary wave as "soliton." 

Solitons in one-dimensional (ID) CA have their own interest and his- 
tory, they have been extensively studied since 1986 by Kennet Steiglitz 
and colleagues as you can see in [56l __TJ [26l [57] . This has been based on 
a variant of classic CA, known as parity rule filter automata (PRFA). 
A PRFA mainly uses newly computed site values as soon as they are 
available and they are analogous to Infinite Impulse Response (IIR) 
digital filters, while conventional CA correspond to Finite Impulse Re- 
sponse (FIR) [56]. Incidentally, yet more sophisticated solitons with 
PRFA were obtained by Siwak in [50] showing large and multiple si- 
multaneous solitonic collisions sequentially (i.e. not parallel mapping). 
We can also see, turbulence solitons in ID CA explored by Aizawa, 
Nishikawa, and Kaneko in [5]. 

Studies of ID soliton CA are important because it allows for fast- 
prototyping of soliton logic. For practical implementations of soliton 
logic, see overview developed by Blair and Wagner in [13] , leads to novel 
designs of optical parallel computers. An interesting implementation 
showing the wave propagation equation in lattice gas simulated with 
a partitioned CA was developed by Margolus, Toffoli, and Vichniac in 
[53] [64], Solitons have found relevant and numerous applications, some 
of them are in: fiber optics, breather waves, non-linear Schrodingcr 
equation, magnets, and recently in proteins and DNA, bio-solitons [T9l 

[Q [51 [sacs]- 

Historically, complex C A have been related to the presence of mobile 
self-localizations (referred as well as: gliders, particles, or waves). The 
most famous CA is the two-dimensional (2D) CA Conway's the Game 
of Life [55], but we can also find a number of samples in ID supporting 
mobile self-localizations, as we can see in [2J [T^l [23 HS1 ISH [SS]. Some 
of them process explicitly signals (not mobile self-localizations) by De- 
lorme and Mozayer in |21) or solve the firing squad synchronization 



-You can find original scanned Russell's papers from Eilbeck's website "John 
Scott Rusell and the solitary wave", http://www.ma.hw.ac.uk/~chris/scott_ 
|russell.html| 1998. 
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problem by Umeo in [SS]. Indeed, we can see how a number of CA 
have been exploited as physical models in [33 EH [68l E] ■ 

This paper is organized as follows. Section [2] gives a general in- 
troduction on CA and basic notation. Section [3] presents experimental 
soliton solutions in CA including solitons in complex ECA rules 54 and 
110, and we will reporting in this paper the soliton reactions emerging 
in Rule 54 from multiple collisions. Later in Sec. 3.3 we displays a 
new ECAM able to solve experimentally the most simple single-soliton 
two-component solution from any initial configuration. Finally we will 
discuss some computing capacities based on solitons. In the last section 
we will discuss final remarks (Sec. [4]). 



| 2. One-dimensional cellular automata 

| 2.1 Elementary cellular automata (ECA) 

A CA is a quadruple (E,<^,/i, cq) evolving on a specific d-dimensional 
lattice, where each cell Xi, i € N, takes a state from a finite alphabet 
E such as x € E. A sequence s € E" of n cell-states represents a string 
or a global configuration c on E. We write a set of finite configurations 
as E™. Cells update their states via an evolution rule ip : E' 1 — > S, such 
that /j, represents a cell neighbourhood that consists of a central cell 
and a number of neighbours connected locally. There are |E| M different 
neighbourhoods and if k = |E| then we have k k different evolution 
rules. 



t + 1 



53 



t + 2 



t + 3 



t + 4 



Figure 1. Dynamic in ECA on an arbitrary one-dimensional array and hypo- 
thetical evolution rule ip. 



An evolution diagram for a CA is represented by a sequence of con- 
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figurations {ct} generated by the global mapping <!> : E n — > E n , where 
a global relation is given as $(c f ) — > c t+1 . Thus cq is the initial config- 
uration. Cell states of a configuration c* are updated simultaneously 
by the evolution rule as: 

ri ■ ■ ■ ) X il • • • ) x i+r) ~ ^ ■ (1) 

where i indicates cell position and r is the radius of neighbourhood y,. 
Thus, the elementary CA class represents a system of order (k = 2, r — 
1) (in Wolfram's notation [35]), the well-known EC A. To represent a 
specific evolution rule we will write the evolution rule in a decimal 
notation, e.g. (finuo- Thus Fig. [I] illustrates how evolution dynamics 
work in one dimension for EC A. 

| 2.2 Elementary cellular automata with memory (ECAM) 

Conventional CA are memory less: the new state of a cell depends on 
the neighbourhood configuration solely at the preceding time step of 
(p. CA with memory are an extension of CA in such a way that every 
cell Xj is allowed to remember its states during some fixed period of its 
evolution. CA with memory have been proposed originally by Alonso- 
Sanz in ! EH EE] . 

Hence we implement a memory function <fr, as follows: 

ct>(x t r\...,x t r\xi)^ St , (2) 

where r < t determines the degree of memory and each cell s, G £ 
is a state function of the series of states of the cell Xi with memory 
backward up to a specific value r. Later to execute the evolution we 
will apply the original rule on the cells s as: 

tp(. . . , Si-!, Si, s j+ i, ...)-> x\ +1 (3) 

to get an evolution with memory. Thus in CA with memory, while the 
mapping ip remains unaltered, historic memory of all past iterations is 
retained by featuring each cell as a summary of its past states from (j). 
We can say that cells canalise memory to the map ip [TO] . 

Let us consider the memory function as a majority memory, 

where in case of a tie given by Ei = Eo from <p hence we take the last 
value Xi. Thus, 4> ma j function represents the classic majority function 
(for three values [50]), then we have that: 

<f) maj {a, b, c) : (a A b) V (6 A c) V (c A a) (4) 
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that represents the cells (x\~ T , . . . , x\) and defines a temporal ring 
s before getting the next global configuration c. Of course, this eval- 
uation can be for any number of values of r. In this way, a number 
of functional memories may be used and not only majority, such as: 
minority, parity, alpha, . . ., etc. (see [TU1 ITTj ). 



CO 



<f>m:r • temporal ring s storing memory 




Figure 2. Dynamics in ECAM on an arbitrary one-dimensional array and 
hypothetical evolution rule ip and memory function <f> m with r = 3. 

Evolution rules representation for ECAM is given in [2H IM1 ED] > as 
follows: 

4>CARm:T (5) 

where CAR is the decimal notation of a particular EC A rule and m 
is the kind of memory used with a specific value of r. This way, for 
example, the majority memory (maj) incorporated in ECA Rule 30 
employing five steps of a cell's history (r = 5) is denoted simply as 
4>R30maj:5- The memory is functional as the CA itself, see schematic 
explanation in Fig. [2j However, computationally a memory function 
has a quadratic complexity calculating its evolution space. 

| 3. Solitons in one-dimensional cellular automata 

A soliton is a solitary wave with non-linear behaviour that preserves its 
form and speed, interacting with some kind of perturbation. The latter 
can be another wave or some obstacle, continuing its travel affecting 
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only its phase and position since each collision. For example a water 
wave travelling and interacting with others waves, they can be found 
also in optics, sound, and molecules [!§] . 

The solitary wave described by Scott become formally represented 
by the Korteweg-de Vries equation [25] as: 

1H + u xxxx + uu x = (6) 

where the function u measure high-wave and x-position at time t, and 
every subindex represent partial differences. Second term represent 
scattering- wave and the last term is the non-linear term |26| . 

However, we will indicate that soliton models related to CA do not 
find some direct relations matching some differential equation solu- 
tions. Nevertheless, Steiglitz has been displayed some properties with 
Manakov systems and PRFA in [32] [57] in the search of computable 
systems collision-based soliton |61j . In addition, Adamatzky in [5J has 
designed a way to manipulate solitons to implement logic gates. On 
the other hand, Chua has developed explicitly an extended analysis 
on how ECA can be described precisely as differential equations and 
cellular complex networks (CCN) in [IB] , 

Although many studies were done in ECA here we cannot find much 
about the soliton phenomena for each rule. Complex ECA are direct 
candidates to explore such reactions from the interaction of their mobile 
self-localizations. Some explorations were described in [6] [8] [12j HH 
IT51 I2TI I5TI IdT] . Solitons in C A are characterized as a set of cells self- 
organized emerging on the evolution space, such complex patterns have 
a form, volume, velocity, phase, period, mass, and shift. Of course, 
not all these mobile self-localizations may work as solitons because 
they depend on its interaction with other structures. Consequently a 
classification is necessary from the evolution space because they cannot 
be inferred from the local rule. 

While a PRFA was designed to yield solitonic collisions calculat- 
ing the new values as soon as they are available, their mobile self- 
localizations present a strong orientation to the left. This is a natural 
consequence of its function to calculate the next cell which evaluate the 
(i — r) f+1 cells |56) . The main and most important difference with con- 
ventional ECA is that those mobile self-localizations working as solitons 
shall be searched explicitly and cannot be deduced to evolve the sys- 
tem. Thus not all complex ECA are able to produce collisions as soli- 
tons although they could evolve some kind of mobile self-localizations. 
Steiglitz has researched amply the PRFA with the goal of reaching un- 
conventional computing devices based-soliton collision, as we can see 

in EmEgnaEa. 

In the next subsections, we will discuss particular cases with com- 
plex ECA and ECAM, displaying exact codifications to get soliton col- 



Complex Systems, volume (year) 1—1 + 



On soliton collisions in complex EC A: Rules 54 and 110 and beyond 



7 



lisions between mobile self-localizations and we will also present some 
computable capacities. In the ECA domain, we have selected and re- 
searched only complex rules 54 and 110, because none other ECA rules 
present a universe with such diversity of mobile self-localizations and 
consequently an ample diversity of collisions as well. In the ECAM do- 
main we will present a single case that solve experimentally the most 
simple singlc-soliton two-component solution from any initial configu- 
ration, the evolution rule 4>R9 ma j-A- 



I 3.1 Solitons in ECA Rule 110 

ECA Rule 110 is a complex cellular automaton evolving with a compli- 
cated system of mobile self-localizations. Its local function is defined 
as follows: 

f 1 if 001,010,011,101,011 
^ no - \ if 000, 100, in ■ (7) 

Figure [3] illustrates the complex dynamics from a typical random 
initial condition selecting the evolution rule (/Jriio- Here we can see 
how a number of mobile self-localizations emerge on its evolution space 
and how a number of them collide. 




See detailed studies of Rule 110 and mobile self- localizations in: 
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glider systenf] 03 S3 EE], universality OH ES 113 El CE3 S3 ED], 
collisions and Rule 110 objects j44j|47]. So generalities can be explored 
from the Rule 110 repository^ 

We will focus on mobile self-localizations that present solitonic re- 
actions. Localizations that have such property are classified in Fig. |4j 
following Cook's notation [T7], here we can see stationary, shift-right, 
and shift-left (displacements) localizations. 

Rule 110 has an unlimited number of collisions as a consequence of 
some extendible mobile self-localizations [JH US] . In this way, first we 
have constructed a set of configurations c coding each localization and 
yielding the solitonic reaction desired. 

To drive collisions and localizations, we will use the set of regular 
expressions 'fi_l' localizations-based to code initial configurations in 
Rule 110, for full details please see [48 |^] 




Figure 4. Set of mobile self-localizations with solitonic properties in Rule 110. 



mobile self-localization 


shift 


period 


speed 


volume 


A 
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2/3 w 0.666666 


6 


B 


2 


4 


-1/2 = -0.5 


8 


Ci 





7 


0/7 = 


9-23 


c 2 





7 


0/7 = 


17 


E 


8 


30 


-4/15 w -0.266666 


21 


F 


4 


36 


-1/9 » -0.111111 


15-29 




14 


42 


-1/3 w -0.333333 


24-38 



Table 1. Mobile self-localizations properties such as solitons. 



Table[T]shows a number of properties to each mobile self-localization 



2 Gliders in Rule 110 http://uncomp.uwe.ac.uk/genaro/rulellO/ 
|glidersRulellO ■ html) 



'Rule 110 repository: http://uncomp.uwe.ac.uk/genaro/RulellO.html 
4 Regular language glider-based in Rule 110: |http: //uncomp . uwe . ac.uk/genaro/ 
rulellO/listPhasesRllO.txt 
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(Fig. |4|, such as: shift, period, speed, and volume (that can be related 
as its mass as well) . All of them shall help us to synchronize collisions 
given a specific phase, where each mobile self-localization may present 
different contact points and collide with others mobile self-localizations. 
To produce a specific collision between mobile self-localizations at a 
given point we must have a full control over initial conditions, including 
distance between gliders and their phases at the moment of collision. 

The notation proposed to codify initial conditions in Rule lfO by 
phases is as follows: 

#l(#2,f ; :-l) (8) 

where #i represents a particular mobile self-localization (given in Ta- 
ble [I]) and #2 represents its phase if it has a period greater than four 
(for full details please see 48 ) . Variable f^ indicates the phase currently 
used where the second subscript j (forming notation U-j) indicates that 
selected master set of regular expressions. 

In |44j we have calculated experimentally the whole set of binary 
collisions between mobile self- localizations in Rule 110, colliding all 1- 
1 mobile self- localizations. Thus in [37J[33] we have reported all soliton 
reactions in Rule 110. 

This way, 18 solitons (between two mobile self-localizations, i.e. bi- 
nary) in Rule 110 can be coded in phases, as follows. 



(a) 


Soliton 1 




yl(fi_l)-6e-G(C,fi_l) — > 


{G,A} 


(b) 


Soliton 2 




Ci(A,f 1 _l)-3e-£(B,fi_l) 


— > {E,d} 


(c) 


Soliton 3 




Ci(A,fi_l)-3e-£(C,fi_l) 


—> {E,d} 


(d) 


Soliton 4 




_F(A,fi_l)-3e-B(f 4 -l) — 


{B,F} 


(e) 


Soliton 5 




C 2 (A,f 1 _l)-3e-£(C,f 1 _l) 


—* {E,C 2 } 


(f) 


Soliton 6 




Ci(A,fi_l)-2e-F(B,f 1 _l) 




(g) 


Soliton 7 




C 2 (A,f 1 _l)-2e-F(A,f 1 _l) 


— > {F,C 2 } 


(h) 


Soliton 8 




A(fxA)-Ae-E(k,hl) — > 


{E,A} 


(i) 


Soliton 9 




A(fi-l)-4e-F(B,fi_l) — > 


{E,A} 


(j) 


Soliton 10 


A(fi_l)-4e-£(G,fi_l) — 


■> {E,A} 


(k) 


Soliton 11 


A(fi_l)-4e-J(H,fiJ.) — 


-> {E,A} 


(1) 


Soliton 12 


F(A,fi_l)-e-£(A,f 1 _l) - 


—> {E,F} 


(m) 


Soliton 13 


F(A,f 1 _l)-e- J B(C,f 1 _l) - 


► {E.F} 


(n) 


Soliton 14 


F(A,fi_l)-e-^(D,fi J.) - 


y {E.F} 


(o) 


Soliton 15 


F(A,f 1 _l)-e-S(E,f 1 _l) - 


-^{E,F} 


(P) 


Soliton 16 


F(G,f 1 _l)-e-£(A,f 1 _l) - 


-^{E,F} 


(q) 


Soliton 17 


F(G,fi_l)-e-F(B,f 1 _l) - 


-^{E.F} 


(r) 


Soliton 18 


F(G,fi_l)-e-F(H,f 1 _l) - 


► {E, F} 
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Figure 5. (a-r) Binary solitons in Rule 110, and one case (s) illustrating 
multiple solitonic collision with six mobile self-localizations, synchronized 
and evolving in 964 generations. 

Of course, from these solitonic binary collisions we can codify and 
synchronize most structures and therefore to get multiple solitonic re- 
actions increasing its complexity. For example, we have the next codi- 
fication: 

(s) Multiple soliton: Ci(B,fi_l)-e-Ci(A,fi_l)-2e-C 2 (A,fi_l)-e- 
F(A,f 1 _l)-e-£(A,f 1 _l)-3e- J B(C,f 2 -l) — ► {E, E, F, C 1 ,C 1 ,C 2 }. 

All these solitons in Rule 110 are displayed in Fig. [5] Each codifica- 
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Figure 6. Pseudo-soliton in Rule 110. 



collisions F 
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F(A,fi-l)-e 


-S(A,fi_l) 


= {A,B,B,F} 


F(A,fi_l)-e- 


B(fi-l) = {B, F} * 


F(A,fi-l)-e 


-B(B,fi.l) 


= {A,2C 3 ,d} 


F(G,fi-l)-e 


-B(fi_l) = {B,F} * 
B(fi-l) = {S 2 ,^ 2 } 


F(A,fi_l)-e 


-B(C,fi_l) 


= {A,C 2 } 
= {C 2 ,A 2 } 
= {A,A\A,E 


F(H,fi_l)-e 


F(G,fi_l)-e 


-S(A,fi.l) 


F(A2)-e-B 


= {B, F} (soliton) 


F(G,fi_l)-e 


-S(B,fi_l) 






F(G,fi_l)-e 


-S(C,fi-l) 


= {B,F} * 






F(H,fi_l)-e 


-S(A,fi_l) 


= {A,C 2 } 






F(H,fi-l)-e 


-B(B,fi.l) 


= {E,A 5 } 






F(H,fi-l)-e 


■B(C,fi_l) 


= {E, A T '} 






F(A2,fi_l)- 


e-B(A,fi_l 


> = {Ci} 






F(A2,fi_l)- 


e-B(B,fi_l 


= {A,B 3 ,E} 







Table 2. Reactions relation between B, B and F mobile slef- localizations in 
Rule 110. 
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tion of (a) to (r) present the binary case, and (s) presents a multiple 
solitonic collision with six localizations, where each is synchronized to 
produce the soliton reaction. 

Yet as an special case in Rule 110, we can find a collision named 
as pseudo-soliton |37j . that works recovering the original localization 
after two collisions. This is performed with B, B and F localizations. 
Localizations B and B have the same period and speed, but its volume 
is different. 

From "Atlas of binary collisions in Rule 110" [H] we have calculated 
the whole set of binary collisions between mobile self-localizations, and 
summarised them in Table [2j We have placed particular attention on 
asterisks labels because they represent precisely the pseudo-soliton in 
Rule 110. Hence we know that reaction F — s- B = {B, F} and also that 
F — » B = {B, F}, thus a loop may be constructed to synchronize such 
collisions. Figure [6] displays such construction, given its codification in 
phases as: 

F(G,f 3 -l)-2 e -F(A,f 1 _l)- e -B(f 1 _l)-5e- J B(B,f 4 _l). 



I 3.2 Solitons in ECA Rule 54 

ECA Rule 54 is a complex CA evolving with an "apparently" sim- 
ple system of mobile self-localizations. Its local function is defined as 
follows: 

/ 1 if 001,010,100,101 , , 

~ { if 000,011,110,111 ' t j 

Figure [7J illustrates the complex dynamics from a typical random 
initial condition selecting the evolution rule fn.54- Here we can see 
how a number of mobile self-localizations emerge on its evolution space 
and collide. Particularity, Rule 54 is able to evolve with emergence of 
glider gun^J patterns since random initial conditions while this fact is 
not common in Rule 110 or from others ECA rules. 

Detailed analysis of various aspects of ECA Rule 54 can be found 
in: localizations systerrj^] [T2l |2"5I 1501 151] , computations [50] , collisions 
[50] , algebraic properties [58J [35] . So generalities can be explored from 
the Rule 54 repository^ 

The set of mobile self-localizations in Rule 54 is significantly small 
compared with Rule 110 (that has a base set of 12 mobile self-localizations). 
Rule 54 has basically four primitive or basic mobile self-localizations 



5 A glider gun is a complex structure that emits periodically a localization, fa- 
mously known in the Game of Life CA. 

6 Gliders in Rule 54 http://uncomp.uwe.ac.uk/genaro/rule54/glidersRule54. 
html 



'Rule 54 repository: http://uncomp.uwe.ac.uk/genaro/Rule54.html 
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O O O OSXCA Evolution 




Figure 7. Random evolution in Rule 54 on a ring of 644 cells to 375 gen- 
erations. White cells represent the state and black cells the state 1, the 
evolution start with a density at 50%. Also a filter is applied to get a better 
view of mobile self-localizations and collisions on its periodic background. 

(stationary, shift-right, and shift- left displacements) and three kinds 
of glider guns [30]. This way, basic mobile self-localizations work to 
produce solitons from multiple collisions. In Fig. [8] we present these 
basic mobile self- localizations following the Boccara's notation |12) . 




Figure 8. Set of mobile self-localizations with solitonic properties in Rule 54. 

Particularly, solitons in Rule 54 cannot emerge from binary colli- 
sions, they are found in multiple collisions. This way, solitons there 
are on the domain of triple collisions and beyond .'!() . Properties for 
these mobile self- localizations are characterized in Table [H 

In this paper, we will report an unexplored set of solitonic collisions 
in Rule 54. They are obtained for systematic analysis by reactions 
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'■£ . 



(b) 



(c) 



(d) 



(e) 



(g) 



(h) 





(I) (m) (n) 

Figure 9. Cataloge of soliton collisions in Rule 54. 



across of multiple collisions. 

This way, Fig. [9] presents 14 kinds of solitons constructed in Rule 54. 
It is easy to recognise that you can derive 28 similar reactions in total, 
because Rule 54 is a symmetric rule and therefore you can obtain the 
next 14 symmetric collisions. 

We can see that Fig. (a) and (b) display two pairs of mobile self- 
localizations producing the same soliton reaction, however the collision 
is different because while in (a) the first pair of mobile self-localizations 
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period 


speed 


volume 
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2 




-2 
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-1 


0-4 


go 





4 





6-2 


9e 





4 





7-3 



Table 3. Basic mobile self-localizations properties in Rule 54. 



delay its trajectory in (b) it advances for six cells. Thus here is possi- 
ble controller intervals of mobile self-localizations trajectories. Similar 
cases are presented in Fig. (c) and (d), but these reactions are between 
three mobile self-localizations. 

Figure (g) starts solitonic reaction with more than four mobile self- 
localizations, but here we employ more space between intervals of mo- 
bile self-localizations (see (f) and (h) as well). Hence we can use several 
mobile self-localizations to preserve the soliton reaction. Noticeably, 
Fig. (i), (j) and (k) display three different kinds of collisions to get 
soliton reactions employing four mobile self-localizations. 

The last set of collisions (Fig. (1), (m), and (n)) display more large 
synchronisations of mobile self-localizations, with different intervals 
and numbers of them. Of course, it is possible to design more so- 
phisticated collisions working with a diversity of packages of mobile 
and stationary self-localizations. 



By the way, recently a soliton is discovered in ECA Rule 26 (Fig. 10 1. 



Figure 10. Soliton in Class II ECA Rule 26. 



I 3.3 Solitons in ECAM 

In this section, we will present the simple single-soliton two-component 
solution [53] for a specific ECAM. The main characteristic is that only 
one mobile self-localization is processed. Thus a mobile self-localization 
with shift-right and shift-left displacement always produce the same 
reaction. 
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Actually an extensive and systematic analysis is done for ECAM, 
in "Designing Complex Dynamics with Memory" [28] . From here, we 
have selected the ECAM rule 0_Rg ma j : 4, because none another rule have 
the same features. 




OSXCA Evolution 



(c) 

Figure 11. Typical snapshots of ECA Rule 9. (a) starts evolution with a 
single cell in state one, (b) presents a random evolution at 50%, and (c) 
other random evolution with small pixels, 360 cells for 331 generations. 

Figure [TT] shows the ECA base that shall be enriched with majority 



memory function (see Sec. 2.2). We study ECA Rule 9 because it 
displays basic interaction of solitons with simple collisions. As we can 
see in (a) a single soliton travels along the evolution space, while in 
(b) we can see a number of interactions during a short history, starting 
from a random initial condition. Extending the evolution space in (c) 
we can observe better how these solitons emerge in ECA Rule 9 and 
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(a) 




Figure 12. Typical snapshots of EC AM rule cj}R9maj-A- (a) starts with an initial 
density of 10%, (b) presents an initial density of 80%, and (c) has an initial 
density of 50%. All evolutions are filtered for the best visualization of mobile 
self-localization interaction, the evolutions are on a ring of 776 cells for 315 
generations. 
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how they collide inside a fast stationary periodic attractor. 

In j29j |38l [39] , we have demonstrated how EC A when enriched with 
memory with memory produces different dynamics. Here we will ex- 
ploit this tool to get simple solitonic reactions. 

Let use the majority memory with r = 4 in ECA Rule 9. Obtaining 
the EC AM rule 4 > R9maj-A, that evolves with two mobile self-localizations 
emerging on its evolution space: S^ fi9moi . 4 = {~0 , *p}- The localization's 
properties are easy to calculate. The ~jf mobile self-localization has a 
volume of 5 x 6 cells, a mass of 12 cells, and moves 2 cells in 5 gen- 
erations (shift-right displacement). While *p mobile self- localization 
has a volume of 5 x 3 cells, a mass of 7 cells, and moves 2 cells in 5 
generations (shift-left displacement). 

Mobile self-localizations emerging in ECAM <pm ma jA preserve the 
solitonic reaction since any collision. Well, here really there are two 
different collisions (two contact points [S5] or phases [25] in every mo- 
bile self-localization) between ~$ and V mobile self-localizations, but 
at the first collision the soliton is preserved because the sequence is 
fused in a string of four cells in state one, while the second reaction 
fuse a string of eight cells in state one. Although finally they open in 
both mobile self-localizations again late of exactly seven generations. 

Thus the automaton 4>Rg ma jA adjusts every string to always evolve 
with the same mobile self-localization and soliton reactions, as follows 
the next relation of collisions: 

~f -> V = and 
t^ < p={ < p,t}. 



Figure 12 illustrates three different random initial conditions where 
the ECAM rule 4>mmayA always evolve in solitonic collisions. First 



evolution (Fig. 12 1) starts with an initial density of 10% for state one, 
the result implies a high production of p mobile self-localizations with 
very few p mobile self-localizations, preserving always the solitonic 
collisions inside bigger fields of p . In the opposite case, the second 



evolution (Fig. 12 d) has an initial density of 80% for state one and 
produce again, high concentrations of V mobile self-localizations with 
some p mobile self-localizations but newly the solitonic reaction is al- 



ways preserved. The final evolution (Fig. 12 :) displays a 50% of states 
one and zero, generating a similar distribution of both mobile self- 
localizations. In all cases, the ECAM rule <f> nama j-A evolve any initial 
condition in solitons. Thus you can begin with any number of mobile 
self-localizations in (j>ssmarA an d the solitons are always produced. A 
characteristic that no conventional ECA have. However, we will men- 
tion that such behaviour can be reproduced identically in other kinds 
of CA, the reversible block CA (or also known as partitioned ID CA 
as well) explored by Wolfram in [69] (chapter 9). 
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3.4 Computing with CA solitons 



Solitons are useful to preserve information such as in the fibre-optic 
communications field. A particular interest is known if such solitons 
could emulate an equivalent Turing machine. Steiglitz et. al have 
designed a number of results trying to reach this goal, please see 
HO |571 ED] |B3], and logic gates with solitons in [5] ITS]. In [35] for 
example, authors have developed a very simple substitution system 
as an implementation of the function addToHead(), based in soliton 
reactions, where also such mechanisms can be designed as a simple 
collider [30]. Figure 13 displays such operations between two mobile 



self- localizations in the EC AM 4>R30maj:8 [291 138j. 

During the last decade we have seen a number of significant advances 
in work with solitons for modelling unconventional computing devices, 
you can see the next references [3], [4] El H31 EH [38], EH EQl [63]. As a 
result, we can see how solitons could be important to develop com- 
putable devices in the construction of equivalent Turing machines. We 
also want to recall the results obtained in ECA Rule 110, where a cyclic 
tag system was developed to perform a computation based-collisions 
with a large number of mobile self-localizations on an incredible global 
synchronization in millions of cells. Solitonic reactions were very useful 
to write binary data and preserve information in the whole mechanism. 
For full details please see [10 H S3 gpj EH] Q 



OSXCA Evolution 




Figure 13. The ECAM , 



>R30maj:8 



presents a solitonic collision that can be coded 



for any n, m £ Z , such that, p% 



QT n „„ „ is always a soliton. 



Figure 14 summarizes the function of solitons in the construction of 



8 Details and large snapshots about cyclic tag system working in Rule 110 |http7] 
//uncomp . uwe . ac . uk/genaro/rulellO/ctsRule 110 .html 
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Figure 14. Solitons working to yield, handle, and control bits in one cyclic 
tag system working in Rule 110. This evolution begins with 793 cells to 
1144 generations, the evolution is filtered suppressing its ether (periodic 
background) . 
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a cyclic tag system in Rule 110. The first important fact is the solitonic 
reaction between packages of A 4 mobile self-localizations versus one E 
mobile self-localization, this collision should introduce bits on the tape 
of the cyclic tag system yielding a series of G% mobile self-localizations, 
also it might also avoid destroying extra E mobile self-localizations that 
cannot infer on the computable system. Besides, when a binary digit is 
into the tape (stationary C% mobile self-localizations) hence a package 
of E mobile self-localizations, coming from the right side shall preserve 
the binary digit and across them with other solitonic collision, because 
forward them self will be transformed in binary data. Here we can see 
how solitons are useful to preserving information and recognising when 
a value will be processed, deleted, or read (for full details please see 

| 4. Conclusions 

We have reported a complete number of solitons in EC A Rule 110 from 
binary collisions. For ECA Rule 54 we have reported a new number 
of collisions that yield solitons, where they could be manipulated to 
develop computable devices, or yet more, complex constructions based 
solitonic reactions. So, finally we have characterized a simple single- 
soliton two-component solution with a simple EC AM (j)R9maj-A, where 
mobile self-localizations always work as solitons even starting from ran- 
dom initial conditions, because each soliton is always constructed from 

<pR9maj-A- 

With regards to memory effect in C A [THl E] . Recently, we have 
studied how a memory function helps to describe dynamics properties 
that are not evident at the first instance [29l GEl |39l [38] . In present 
paper, the majority memory selected in ECA Rule 9 opens a new evo- 
lution rule ECAM 4>Rg ma j-A able to simulate solitons. Of course, frag- 
ments of the original evolution rule determine such dynamics. 
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